Field-induced XY behavior in the S = ^ antiferromagnet on the square lattice 



Alessandro Cuccoli/'^ Tommaso Roscilde/'^ Ruggero Vaia,'^'^ and Paola Verrucchi^^ ^ 

^ Dipartimento di Fisica dell'Universita di Firenze - via G. Sansone 1, 1-50019 Sesto Fiorentino (FI), Italy 
^Istituto Nazionale per la Fisica della Materia - U.d.R. Firenze - via G. Sansone 1, 1-50019 Sesto Fiorentino (FI), Italy 
^Istituto di Fisica Applicata 'Nello Carrara' del Consiglio Nazionale delle Ricerche - via Panciatichi 56/30, 1-50127 Firenze, Italy 

(Dated: February 2, 2008) 

Making use of the quantum Monte Carlo method based on the worm algorithm, we study the 
thermodynamic behavior of the S — ^ isotropic Heisenberg antiferromagnet on the square lattice in 
a uniform magnetic field varying from very small values up to the saturation value. The field is 
found to induce a Berezinskii-Kosterlitz-Thouless transition at a finite temperature, above which 
a genuine XY behavior in an extended temperature range is observed. The phase diagram of the 
system is drawn, and the thermodynamic behavior of the specific heat and of the uniform and 
staggered magnetization is discussed in sight of an experimental investigation of the field-induced 
XY behavior. 
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Field-induced effects in low-dimensional antiferromag- 
nets have been the subject of renewed interest in the 
last few years; while on the experimental side fields of 
very high intensity have become available, on the the- 
oretical side the possibility of inducing novel magnetic 
phases via application of a strong field have been pointed 
outfliii. 

In this paper we consider the two-dimensional quan- 
tum Heisenberg antiferromagnet (2D QHAF) in a uni- 
form magnetic field, described by the Hamiltonian 

i,d i 

where i — (ii, 12) runs over the sites of a square lattice, 
d connects each site to its four nearest neighbors, J > 
is the antiferromagnetic exchange coupling, H is the ap- 
plied Zeeman field, and [S'p = S'(5'-f 1). We will hereafter 
use reduced temperature and magnetic field, t = T/ J and 
h = gn^H/{JS). 

The rich phenomenology l5l| of the model is ruled by the 
interplay between the exchange and the Zeeman terms 
in Eq. (0). The apphed field breaks the 0(3) symmetry 
of the isotropic model and induces a uniform alignment 
in the z direction; such alignment frustrates the antifer- 
romagnetic order along z but does not clash with an- 
tialignment on the xy plane, where 0{2) symmetry stays 
untouched. For infinitesimally small fields one hence ex- 
pects the spins to lay antialigned on the xy plane, and 
progressively cant out of it as h is increased. Satura- 
tion occurs at the critical value he = 8, above which the 
ground state displays uniform ferromagnetic alignment 
along the z direction. In the range <h<hc one may 
also expect thermal fluctuations of the z spin components 
to be smaller for larger field, while no such a reduction 
should occur as far as the x and y components are con- 
cerned. The above picture clearly suggests the model to 
share essential features with the easy-plane 2D QHAF Q . 

In the classical limit {S-^oo, JS^^Jc\, 



g^^H ^ h Jc\) both analytical 0, and numerical 
calculations revealed the occurrence of a Berezinskii- 
Kosterlitz-Thouless (BKT) transition [H [13, for all 
values of h below saturation. In the quantum S = ^ case, 
evidence of a field induced BKT transition was recently 
achieved for small fields by means of quantum Monte 
Carlo (QMC) simulations based on the continuous-time 
loop algorithm Unfortunately the loop algorithm 
loses its efficiency exponentially as the field and/or the 
inverse temperature are increased [l^ : this fact has 
so far prevented from a systematic investigation of the 
strong-field regime |l3|. However, the recently proposed 
worm Q (or directed-loop |l4|) algorithm , working in 
any lattice dimension, is effective also in the presence of 
a uniform field of arbitrary intensity. Our QMC simula- 
tions are in fact based on such an algorithm, which is a 
pure-quantum cluster algorithm (by this meaning that 
it has no classical analogue, so far), that takes the field 
into account by the dynamical nature of the process 
of cluster (worm) growth, with the possibility that the 
worm traces back its route after bouncing at some point. 
In particular the update process is irreversible (i.e., the 
inverse of a single update step with finite probability 
can have a vanishing probability), thus reflecting the 
time-reversal symmetry breaking due to the presence 
of the field. A detailed description of the algorithm is 
beyond the scope of this paper and can be found in 
Ref. Hsl Our simulations were performed on a L x L 
square lattice (L = 16, 32, 64, and 96), each consisting 
of 10* MC steps for thermahzation and of (1-1.5) xlO^ 
MC steps for evaluation of thermodynamic observables. 
During thermahzation, the number of worms produced 
at each step is adjusted so that the total length of the 
worms in the imaginary-time direction roughly equals 
the size of the (2-|-l)-dimensional lattice, L'^/t; this 
number is then kept fixed during the measurement 
phase. In this way, autocorrelation times of the order of 
unity are achieved for all values of the field. 
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TABLE L t^KT W ^ obtained by finite-size scaling analysis. 
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0.1 


0.175(5) 
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0.254(5) 


4.0 


0.282(5) 


0.2 


0.195(5) 


2.0 


0.292(5) 


6.0 


0.202(5) 


0.4 


0.213(5) 











Thanks to this very effective tool, we got access to the 
thermodynamic behavior of the model (Q) with S — ^ 
and h varying from to he, looking for signatures of 
field-induced BKT behavior that can be object of exper- 
imental observation. In particular, we have focused our 
attention on the specific heat c{t, h) and the field-induced 
uniform magnetization m^{t,h), which are easily acces- 
sible to experiments, as well as on the staggered mag- 
netization along the field axis m^j^{t,h), which provides 
further insight in the microscopic ordering mechanism. 
Seven values of the field have been considered: h = 0.1, 
0.2, 0.4, 1, 2, 4, and 6. 

First of all we have performed a detailed finite-size 
scaling analysis in order to check whether or not the pre- 
dicted critical scaling behavior of the in-plane staggered 
susceptibility and of the helicity modulus are re- 
produced at consistent temperatures: having got a posi- 
tive answer, we may state that a finite-temperature phase 
transition of BKT type occurs in the model for all the 
considered field values. The estimate of the critical tem- 
perature corresponding to each field has been obtained 
via the same procedure used in Ref. for the easy-plane 
2D QHAF: the resulting values are listed in Tabled they 
are consistent with previous results Q for h = 0.2 and 
h = OA. In Fig. we report the phase-diagram of the 
model for S = and S = oo, the latter as from classi- 
cal MC results 1^. We observe that the effect of quan- 
tum fluctuations is limited to a strong renormalization 
of tg^rj. with respect to the classical case, but the field 
dependence is qualitatively the same. 

Let us first concentrate on the weak-field regime: 
from various analytical arguments |^ ItR one can 
infer a mapping of the 2D QHAF in a field on a 
weakly easy-plane magnet with exchange anisotropy, i.e., 
Hamiltonian |^ with H ^ and the additional term 
— J'^J2id'^i'^i+d'^ effective exchange anisotropy A 
is expected to scale with the field as A oc /i^. Hence, the 
BKT critical temperature is expected to obey, both in 
the classical and in the quantum case 0|, the following 
expression 



Ah) 



ln(C//i2) 



(2) 



where is the spin stiffness of the isotropic model, and 
C a constant. Fitting our results for the lowest three val- 
ues of h to the above expression we get ~ 0.22 J, and 
excluding the third value we find ~ 0.19. Remarkably, 




FIG. 1: Phase diagram of the S — ^ 2T) QHAF in a magnetic 
field. Open symbols refer to the classical limit of the model, 
from Ref. 13 the triangle is a QMC resuh from Ref.[l3 Inset: 
*BKT (h) and tl (h) (see text) for weak fields; solid and dashed 
lines are logarithmic fits to the first three points of each data 
set, and the shaded area marks the region of disordered XY 
behavior. 



these values are consistent with the renormalized spin 
stiffness of the S'=i isotropic 2D QHAF, =0.180 J 
[igj . We also notice that the fitting curve keeps interpo- 
lating the data up to h^2. Finally, by directly compar- 
ing the low-field behavior of t^^,^ (h) with that of t^,^^ (A) 
in the weakly planar antiferromagnet 0, we obtain an 
excellent agreement for A « (3/8) /i^. 

Moving towards higher fields two different effects are 
expected: (i) the fluctuations of the z components be- 
come smaller, resulting in an enhanced effective easy- 
plane anisotropy, and (ii) the average projection of the 
spins in the xy plane decreases, due to the increasing 
uniform magnetization. Globally, the system behaves 
as a renormalized planar rotator with progressively re- 
duced rotator length. It is remarkable that, despite the 
spin configurations being characterized by a smaller and 
smaller projection in the xy plane, the XY character is 
apparent even for h close to he and the transition ver- 
ifies all the predictions of the BKT theory. The inter- 
play between the two competing field effects yields the 
non-monotonous dependence of tg^rj,{h): for low field the 
reduction of z fluctuations is dominant and ^b^t (^) 
creases with h, starting as in Eq. while for higher 
field spin canting prevails and tg^^{h) decreases, even- 
tually vanishing at the saturation field he- Therefore, a 
maximum in ig^T (^) connects the two limiting behaviors, 
as already observed in the classical phase diagram: the 
comparison in Fig. ^ shows that for S=^ the maximum 
shifts to slightly higher field as a consequence of quantum 
fluctuations. 

We now consider the temperature dependence of some 
relevant observables, beginning with the specific heat 
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FIG. 2: Ac{t,h)/t = [c{t,h) - c{t,0)]/t vs temperature, 
for four different field values. Insets: magnetic specific heat 
c{t,h) compared to c{t,0) (thick solid line). The zero-field 
specific heat c{t, 0) is obtained by interpolating numerical and 
analytical data from Refs. [g^, |M [g^); the arrows mark the 
estimated BKT critical temperatures. 

c{t,h). While our data resolution for /i<0.4 prevents 
from observing significant deviations with respect to 
the zero-field system, the results for the four largest 
fields are shown in Fig. |21 in particular, what is plot- 
ted is the specific heat variation upon application of the 
field, Ac{t,h) = [c(t, /i) - c(i,0)], divided by t. This 
quantity equals the difference of the entropy derivatives 
dtS{t,h) — dtS{t,0) and allows us to draw the follow- 
ing picture: At low temperature the entropy increase 
is smaller than in zero field reflecting the presence of 
quasi long-range order induced by the field via stabiliza- 
tion of bound vortex/antivortex pairs. Slightly above the 
BKT transition {t about 20 30 % larger) a sharp en- 
tropy increase occurs, which we interprete as due to vor- 
tex unbinding. When the temperature is further raised 
all vortices are free and then S{t, h) increases slower than 
0); eventually the entropy difference vanishes in the 
fully disordered system at t oo. For incresing field the 
peak of c{t, h), which mimics the BKT peak of the XY 
model, moves to lower temperature, thereby getting nar- 
rower, as shown in the insets of Fig. [3 It is worth noting 
that Ac{t,h), which bears very clear signatures of the 
BKT behavior, is also easily accessible to experiments, 
because non-magnetic contributions cancel in its defini- 
tion. 

We have also considered the field-induced uniform 
magnetization, m^{t,h) = (S'|); such a quantity, which 
can be experimentally determined via standard magne- 
tometry measurements, is also a highly precise output of 
the QMC simulations. In Fig. Owe report m^{t, h)/h for 
different fields. For h^2 and for high enough tempera- 
ture this quantity is found to coincide with the uniform 
susceptibility of the zero-field case, thus showing that 
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FIG. 3: Field-induced uniform magnetization m^{t,h) vs 
temperature, for different field values. The stars represent 
the zero-field uniform susceptibility from Refs. l2Ctl22l: arrows 
as in Fig. |5| 
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FIG. 4: Finite-size staggered magnetization mlj^{t,h) vs 
temperature, for different field values. The x's represent the 
same quantity for the S =^ 2D XY model. Arrows as in 
Fig. El 

the magnetization process is linear. Upon lowering t the 
nonlinearities show up and the uniform magnetization 
changes completely its temperature dependence, display- 
ing a minimum at t = t^{h). This feature marks the onset 
of XY behavior: as the temperature is lowered below t* , 
the system is increasingly magnetized along the z-axis 
and the short-range antiferromagnetic correlations of the 
z spin components, as well as their thermal fluctuations, 
are suppressed, in turn stabilizing the canted configura- 
tions. A remarkable feature is that this crossover to XY 
behavior is located at a temperature f^ih) well above the 
critical point. For ft. ^2 the minimum above ig^T disap- 
pears, and the most prominent feature is rather the shift 
of the broad maximum to lower temperature. 

The crossover from isotropic to XY behavior in low 
fields {h < 2) can be also detected in the temperature de- 
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pendence of tlie finite-size staggered magnetization along 
the hard z-axis HI, m|^(t,/i) = Ei(-)'5'||), 
shown in Fig. 0] In absence of long-range order mf r is 
known to scale to zero as 1/L (for large enough L) [2J|, 
and j^L is hence a bulk property of the system. In the 
limit t (X the system behaves as a collection of para- 
magnetic spins so that, by applying the central-limit the- 
orem, one finds L ^/2S{S + l)/(37r) « 0.399 for 
S — ^. As t decreases, ml L in the isotropic 2D QHAF 
monotonically increases and diverges for t-^i). In the 
2D XY model the same quantity decreases below the 
infinite-t value, due to the suppression of out-of-plane 
fluctuations. The coexistence of both the above behav- 
iors is most clearly observed for and the appear- 
ance of a maximum in ^ at t = t* also marks the 
crossover from isotropic to XY behavior. 

Hitherto, we have identified two ways of locating a 
crossover temperature, at least for low and interme- 
diate fields. One can check that for low fields these 
temperatures match with each other, so 
that the crossover temperature is unambiguous and 
its estimates are i*(0.1) = 0.22(1), t*(0.2) =0.25(1) and 
t*(0. 4) = 0.29(1). In this regime the crossover tempera- 
ture is expected to follow a logarithmic behavior analo- 
gous to that of Eq. lU (with a different coefficient C") 
as suggested in Refs. 6, and,23i, a fit of t*{h) for /i<0.4 
gives /9c. = 0.19, again in good agreement with the known 
value 'l3| of the spin stiffness for S=-^. For larger 
fields t*{h) is systematically higher than i*(/i), suggest- 
ing that the crossover phenomenon extends over a wider 
temperature range, but we note that for intermediate 
fields the same fitting function still interpolates the data, 
t*(l) = 0.36(1) and t*{2) = 0.48(2), as shown in the inset 
of Fig. Finally, for strong fields the explicit signature 
of the crossover gradually disappears from * L, and 
for h = 6 this quantity is nearly monotonic as in the XY 
model; the antiferromagnetic interaction of the z compo- 
nents is then almost completely overcome by the applied 
field. 

In conclusion, we have studied the S — ^ two-di- 
mensional quantum Heisenberg antiferromagnet on the 
square lattice in an arbitrary uniform field by means of 
the quantum Monte Carlo method based on the worm al- 
gorithm. Our results point out that an arbitrarily small 
field is able to induce a BKT transition and an extended 
XY phase above it, as in the case of an easy-plane ex- 
change anisotropy. The field-induced XY behavior be- 
comes more and more marked for increasing fields, while 
for strong fields the antiferromagnetic behavior along the 
field axis is nearly washed out, so that the system be- 
haves as a planar rotator model with antiferromagnetism 
surviving in the orthogonal plane only; the BKT critical 
temperature vanishes as the field reaches the saturation 
value he and the effective rotator length goes to zero. The 
model in a moderately strong field represents an ideal 
realization of the XY model: XY behavior can be de- 



tected by measuring standard non-critical quantities, as 
the specific heat or the induced magnetization; this opens 
the possibility for an experimental realization of the XY 
model in purely magnetic systems, and for a systematic 
investigation of the dynamics of vortex/antivortex exci- 
tations. 
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